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Abstract 

We have developed a multi-phase SPH method to simulate arbi- 
trary interfaces containing surface active agents (surfactants) that lo- 
cally change the properties of the interface, such the surface tension 
coefficient pQ. Our method incorporates the effects of surface diffusion, 
transport of surfactant from/to the bulk phase to/from the interface 
and diffusion in the bulk phase. Neglecting transport mechanisms, we 
use this method to study the impact of insoluble surfactants on drop 
deformation and breakup in simple shear flow and present the results 
in a fluid dynamics video. 

1 Introduction 

Exposing drops to extensional flows such as e.g. a simple Couette flow, the 
viscous forces along the interface tend to deform the drop and elongate it to 
an ellipsoid-type shape. The balancing force to stop the deformation due to 
the shearing is the surface tension. When two pure fluids of different types 
are in contact, the resulting surface tension force is only proportional to the 
local curvature and is normal to the interface. Depending on the strength 
of this force and the viscosity ratio between the two fluids A, drops are 
deformed to a steady ellipsoid shape or break up. The correlation between 
the breakup behaviour and the flow parameter expressed by the Capillar 
number Ca (A) is known as the Grace curve [6]. 

Adding surface active agents (surfactants) to a multiphase system can 
strongly alter the flow phenomena. Surfactants mainly affect the surface 



tension coefficients between two fluids when replacing fluid molecules at the 
interface with surfactant molecules. Hence, surface tension gradients along 
the interface can occur, resulting in Marangoni forces [8]. 

Here, we only focus on the case of insoluble surfactants, i.e. surfactants 
are initially added to the interface and cannot dissolve to the adjacent fluid 
phases. Bazhlekov et al. [2J studied the effect of insoluble surfactants on 
drop deformation and breakup in simple shear flow with a boundary-integral 
method and clearly describe the different breakup modes. However, due to 
the nature of their method, an interface capturing scheme is required and 
breakup detection requires special procedures. By the use of a Lagrangian 
particle method we can avoid these steps and can handle strong interface 
deformations naturally. 



2 Governing equations 

Following the well-known weakly compressible SPH method [7], we use an 
equation of state (EOS) to relate the pressure to the density and solve the 
isothermal Navier-Stokes equations on a moving Lagrangian frame. Formu- 
lating the surface force as the gradient of the surface stress tensor [3] 

= V • [a (I - n ® n) <5 S ] = - (am + V s a) <fe , (1) 

surface tension can be split into a normal Capillary force aKn<5s with the 
curvature the normal vector of the interface n and the surface delta 
function 5% and the tangential Marangoni force V s aS^ (V s is the surface 
gradient operator V s = (I — n (g) n) V). 

Assuming insolubility in the phases, the evolution of surfactant on the 
interface is governed by a diffusion equation 

^ = v s • D s v s r , (2) 

at 

where T and D s are the interfacial surfactant concentration and the diffusion 
coefficient matrix (in case of isotropic diffusion D s = D s ■ I), respectively. 

To close our model, we relate the interfacial surfactant concentration 
r to the surface-tension coefficient a by a constitutive equation. Widely 
used in literature, the Frumkin isotherm or the Langmuir model [9] are 
known to agree reasonably well with experimental data. By the purpose of 
visualization, however, a simple linear relation between a and T is sufficient 
and is employed here. 
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3 Numerical results 



The results of our simulations are presented in |Video~T| 

As a reference, in the first part of the movie we show the effect of the 
Capillary number on the droplet behaviour in the shear flow. At sub-critical 
conditions, the drop deforms to a steady ellipsoid. When the shear forces 
are strong enough, the surface forces cannot counterbalance the shearing of 
the drop and breakup occurs. 

When surfactants are present at the interface, another very complex phe- 
nomenon can be observed, namely the so-called tipstreaming. At the tips of 
the deformed drop the interface is compressed and consequently, the surfac- 
tant concentration increases locally. Thus, surface tension gradients develop 
and Marangoni-forces change the behaviour of the drop dramatically. As 
the surface tension is very low at the tips of the drop, very sharp tips with 
a high curvature develop. Under critical conditions, this effect results in a 
singularity where the tipstreaming starts. In our simulations we can show 
that the thickness of the stream is of the order of the size of a particle. 

At the end, we show the possibilities to influence the tipstreaming by 
changing the strength of the surface diffusion. When diffusion is strong, 
concentration gradients are smoothed and the Marangoni-forces are weaker. 
Hence, the curvature at the tips can be reduced and tipstreaming is sup- 
pressed. Another way of controlling the tipstreaming is by changing the 
properties of the surfactant, which is not shown here. 
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